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Abstract 



It is shown that a SU(1,1) algebra may be used to provide a unified description of the 
simple hamonic oscillator and the angular momentum algebras and a class of other semi- 
infinite algebras. A normal ordered representation of a Unitary operator U constructed from 
the generators of a SU(1,1) algebra, which is a generalisation of the Baker- Campbell - Haus- 
dorff relation for Lie algebras, is given. It is shown that the normal ordered representatiion 
of U may be used to calculate expectation values which are functions of the parameters used 
to construct the operator. The functions so constructed satisfy certain recurrence relations 
and the entire set of functions may be interpreted in terms of diagrams similar to the Pascal 
triangle for binomial coefficients. Coherent states, squeezed states and rotation matrices of 
the angular momentum algebra emerge as special cases. 



1 Introduction 

We consider a set of operators defined on a set of orthogonal basis vectors of a separable, finite or 
infinite dimensional complex vector space H with vectors \j), vifhere j = refers to the vacuum 
vector, positive integer values of j refer to 'particle' vectors and negative integer values of j refer 
to 'hole' vectors. Let the matrix representations of the operators L, its hermitian adjoint = R 
and their commutator S = [L, R] in the vector space H be defined by 



Ljk = Xj Sj^k-^i , Rjk = \k Sj^k+i , j, A; = 0, ±1, ±2, .... (1) 



- (A^ -X^_,) 5,k (2) 



where Xj are real numbers. The operators i?, L and S are the raising, lowering and number 
operators of the algebra. Since 

[L,S],k - A, (A^+i + A2_i - 2A2) 5,^u-i (3) 

the algebra closes if the parameters Xj have j dependence of the form 

X) - cj {a + j) {p + j) (4) 

so that 

X] - X]_^ = a{2j - I +a X)+^ + X)_^ ~ 2 X] = 2 a (5) 
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where a, /? and a are real numbers. The resulting closed algebra 

\R,L\ ^ S , \L,S\ = 1 a L , \S,B\ = 2 a R (6) 
can be identified as an SU(1,1) algebra. 

Of the three parameters present in the algebra cr is a scale parameter while a and /3 provide 
different possibilities for the realisation of the algebra depending on whether (a, /?) are positive 
or negative, integer or non-integer. It is possible to realise closed or infinite or semi-infinite 
algebras by suitable choices of the parameters. Three familiar realisations of the algebra in © 
can be identified. 

• The creation and annihilation operators and a of the Simple Harmonic Oscillator ( hereafter 
refered to as SHO) obeying the commutation relation [a, — 1 when 2i — and 25* = (a^a -|- 
aa^) is an example of the semi-infinite realisation of the algebra. This possibilty corresponds to 
choosing either a = 2/3 = 1 or 2a = /3 = 1 and either all the even number particle states or all the 
odd number particle states of the SHO provide an example of a realisation of the algebra. 

• Ifi=:a, R — a\ a = l,/3^-oo,cr— T^O and (t/3 ^ 1 the resulting algebra is again the 
algebra of the SHO and all the number states of the oscillator taken together provide another 
semi-infinite dimensional realisation of the algebra. 

• A finite dimensional example is the angular momentum algebra whose generators \plJ± = J^; ± 
iJy and Jz satisfy the commutation relation \3x^ Jy\ = iJz and its cyclic permutations and can 
be mapped on to the RLS scheme by R = J+, L = J^, S — Jz with a — a = n, (3 ~ —n~l 
and different choices of integers n provide different finite dimensional realisation of the algebra 
corresponding to different values of = n(n -|- 1). Half- integral angular momenta can also be 
accomodated in this scheme by switching to the basis vectors |iV) with N taking only half- integral 
values instead of integral values and finite dimensional realisations are possible now whenever 
a = TV and ^ = -iV - 1. 

In this study we consider general cases of the realisation of the SU {1,1) algebra when q;,/3 
and a take arbitrary values. When a and f5 take values which are different from integers or 
half-integers infinite dimensional algebras arise and these algebras may prove useful in new 
contexts. The SU(1,1) algebra has played an important role in Mathematical Physics. Various 
representations of SU(1,1) have been extensively studied (Barut and Fronsdal 1965, Holman and 
Biedenharn 1966). The SU(1,1) group and its Schwinger representation (Schwinger 1965) have 
been used in Quantum Optics in connection with the study of parametric amplifiers (Louisell 
1977), interferometers and squeezed states (Yurke et al 1986). 

2 Normal Ordered Representation 

It is well known in Quantum Optics that the squeezed states and coherent states of light may 
be generated from the vacuum state of a simple harmonic oscillator mode by the action of a 
Unitary operator acting on the vacuum state. The generators (i?, L, S) may be used to construct 
a Unitary operator Ui — exp{iy(R + L)) which is a generalisation of the Unitary operator 
arising in the study of the generation of squeezed states and coherent states of light from the 
vacuum. The normal ordered representation, in which the lowering operators are aligned on the 
right and the raising operators are aligned on the left and the anti-normal ordered representation, 
in which the raising operators are on the right and the lowering operators are on the left, of Ui 
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may be found using the method discussed by Truax (1985) and Sukumar (1989) and given in the 
form 

C/i = (exp(z2;(i? + L))) = (exp(z2//i?)) {g^^) (exp(z2//L)) (7) 

3 

= {eMwfL)) (5](|j) {g-^^) (j|) (exp(*y/i?)) (8) 

3 

tanhw-i/CT , ^ , „N 

/ = p— , .9 = sech y^/^ , p, = {2j-l + a + /3 . 9 

This result is a generahsation of earher resuhs which were derived for special values of the 
parametrs a, /3, a and is generalisation of the Baker- Campbell Hausdorff relation. The set of 
functions Gn defined by 

G„ = (-^)" {n\ C/i |0) , n = 0,l,2,... (10) 
may be evaluated using the anti- normal ordered representation of C/i to get 

G„(a,/3;a;y) = (-z)" £ ^^^^^ (nl^^^'-^ |j) g^''^=~'^~^'> (j|i?^|0) ^ . (11) 

j=n 

The matrix representations of the generators {R, L) and the power series expansion of the hy- 
pergeometric functions of the second kind (Gradshteyn and Ryzhik 1965) may be used to show 
that 

G„(a,/3;a;2/) =A„(^^^^i^^)"(sechyV^)"+^-i 2^1(1 - a, 1 - /3; 1 + n; ~ sinh^ y^^) (12) 

An 



' " Y r(a)r(7i+l) r{l3)T{n + l) n! f}. 

Using L\n + 1) = A„|n) and R\n) — A„|n + 1) and (|10D it can be established that G„ satisfy the 
recursion relation 

d 

-T^Gn+l ~ (A„ Gn — A„+i G„+2) ■ (14) 

This recursion relation may also be proved directly using the properties of the hypergeometric 
functions. It is evident that the G„ are symmetric in a and j3. It is clear from (jl2p that a 
variety of solutions are possible depending upon the values of a, /? and a. We note that the 
functions G„(j/) are the expansion coefficients that describe a superposition of states arising from 
the action of the Unitary operator Ui acting on the vacuum state and they depend upon the 
generators of the algebra and the parameter y. 



2.1 Diagrammatic Representations 

The power series representation of the set of solutions Gn may be displayed using a set of nodes 
on a grid set in the form of a rhombus. Such diagrams in general correspond to generating a 2- 
dimensional array of points (j, k) such that the sites for which j+k is an even integer are occupied. 
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Weight factors Ajk can be associated with the nodes such that Ajk = ii j + k = 2m + 1. Using 
such an array it is possible to study a class of functions of two variables: 

oo j oo CKJ OO OO 

j=0 k=0 k=0 j=k j=0 k=0 

where the sum over the functions fj corresponds to constructing functions based on the weights 
associated with the nodes in a row and then summing over such functions with weights 
while the sum over the functions fj corresponds to constructing functions based on the weights 
associated with the nodes in a column and summing over such functions with weig hts g'=. The 
nodes in the diagrams then code for a particular class of functions of two variables such that the 
rows define a sequence of polynomials and the columns define a sequence of functions with infinite 
number of terms in their power series. Diagrams of this kind are relevant for the algebra under 
discussion for the following reason. The matrix representations of {R, L) considered in this paper 
have non- vanishing elements only on either the diagonal above or the diagonal below the principal 
diagonal. Consequently the k^^ powers of these matrices have non- vanishing elements only along 
the fc*'* diagonal above or below the principal diagonal. Ui contains all powers of the matrix 
M = R + L and the matrix elements of powers of M naturally lend themselves to representation 
in terms of diamond like formations. When the 'hole' states decouple from the 'particle' states 
one half of the diamond formation is adequate and leads to triangular diagrams. We examinbe 
this case first. We consider the case when one of (a, /3) equals 1 and the other takes the positive 
value p to identify some of the features of the scheme introduced in this section. 



2.2 Study of = (l,p) , a = 1 

For this case A| = (1 + j){p + j) which means that A_i = and as a consequence the matrix 
representations of {R, L) exhibit a complete decoupling of the 'particle' and 'hole' states and it 
would be adequate to restrict attention to the set of states \n) with n in the range [0, oo]. The 
hypergeometric function in (|12l) simplifies to just 1. In this case a power series representation of 
the solutions (IT^ in the form 

G.iv) - E jtr^. (16) 

m=0 ^ ' 

leads to the recursion relation 

Cn+l,m — A„ Cn,m + ^n+l Cn+2,m-l (17) 
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which may be interpreted in terms of a diagram of the form 



Ao 








• 

Ai 

\ 




Ao 
\ 


Ai 


• 

A2 
\ 


Ao 


• 

Ai 

\ 


A2 


Ao 

\ 


Ai 


• 


Ao 


• 





where the nodes in the columns represent Cn.m for fixed n and varying m. The vahies attached 
to the nodes correspond to the sums arising from the addition of aU contributions from paths 
pointing towards the node, with each path contributing an amount equal to the product of the 
value of the node at the start of the arrow and the weight factor attached to the arrow. The 
power series expansion of 



/ T(n + v) 

G„(l,p;l;t/) = Gn{p,l;l;y) = \ L, V^,^; , scch^t/ tanh"t/, n = 0,1,2,... (18) 

is related to the values at the nodes of the columns of the diagram given above. The presence of 
square root factors in these functions may be removed by defining 



'r(p)r(n 


+ 1) 


r(n + 


p) 


' T{n + 


p) 


r(p)r(n 


+ 1) 



Gn{p;y) = J^^^TT^^^ Gr.il, p;l;y) = sech^' y tanh" y (19) 

T{n+p) 
T{p)T{n + l) 



Gn{p;y) = Gn{l,p;l;y) = Jt.^} sech^ y tanh"y (20) 



which satisfy 

— Gn+i{p;y) = (n + l) Gn{p;y) - (n + 1 +p) G„+2(p;y) (21) 

■^Gn+i{p;y) = (n + p) Gn{p;y) - (n + 2) G„+2(p;y) (22) 
oy 

We note that the case p = 1/2 leads to Squeezed states in Quantum Optics and in particular 
n = corresponds to the squeezed vacuum in Laser Physics (Yuen 1976). The scheme described 
above can be illustrated with explicit evaluation of the values attched to the nodes for specific 
values of p. 

For example if p=l for which \j = (j + 1) and 

G„(l;y) = G„(l,l;l;y) = sech y (tanhy)" , n = 0,l,2, (23) 
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The expansion of these functions in powers of y leads to integer coefficients after factorial factors 
are extracted and the set of integer coefficients may be captured in a scheme of the form 

1 



1 



5 



1 

\ 










1 






1 




2 








\ 










2 


1 




2 


3 


\ 






\ 




5 






1 




2 


3 






\ 










28 


1 




2 




\ 










61 






1 

















6 



61 

which is reminiscent of the Pascal triangle for the binomial coefficients. The weight factors 
associated with the links are simply the sequence of XjS when viewed along the diagonals either 



going down the right or down the left. In this diagram the numbers at the nodes of the different 
columns relate to the power series representations of G„ as given by 

Go(l;y) = sechy = 1- ^+ 5^-61^ + .. (24) 

3 5 

Gi{l;y) = sechj/ tanhy = y - 5 |- + 61 |- - .. (25) 

G2(l;y) = sechy tanh^y = 2|- - 28 ||- + .. (26) 

In general the n*'* column of the triangle scheme can be related to the power series for G„_i (1; y). 
It is simple to verify that 

^G„+i(l;y) = (n+1) G„(l;t/) - (n + 2) G„+2(l; y) . (27) 

The numbers at the nodes in the first column of the diagram are the Euler numbers. To find 
higher order terms in the power scries expansions the diagram can be extended by proceeding 
further along the outer diagonal thereby generating a larger diagarm. 

As the next simple we example consider the case when p=2 for which Xj = {j + l)(j + 2). Now 
the resulting functions 

G„(l,2;l;y) = G„(2,l;l;y) = V^H^: sech^ y (tanhy)" , n = 0,l,2, (28) 

satisfy the recursion relation 

^G„+i = Vn + 2 (Vn+1 Gn - Vn + 3 Gn+2) ■ (29) 
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The factor vn+T present in these functions may be removed if we define 

Gn{2;y) = -yl= C?„(l,2;y) = sech^ y tanh"?/ (30) 
Vn + 1 

Gn{2;y) = x/^H^: G„(l, 2; y) = (n + 1) sech^ y tanh" y (31) 
which satisfy the recursion relations 

^ G„+i(2 : y) = (n + 1) G„(2 : y) - (n + 3) G„+2(2 : y) (32) 

1^ G„+i(2 : y) = (n + 2) (g„(2 : y) - G„+2(2 : y)) . (33) 

The expansion of G„ and G„ in powers of y lead to integer coefficients after factorial factors are 
extracted and the resulting set of integer coefficients may be characterised by triangular schemes 
of the form 

1 1 



16 



1 

\ 








2 

\ 










1 








2 






2 




2 




1 




3 








\ 








\ 










2 


2 






6 


1 




3 


3 


2 




2 


4 


\ 






\ 


\ 






\ 




8 




6 




16 






2 




2 


4 


1 




3 


3 






\ 








\ 










40 


16 






120 


1 

\ 




3 




2 

\ 




2 






136 








272 






2 








1 

























272 272 

The numbers at the nodes of the different columns of the diagram on the left relate to the power 



series representations of G„(2 : y) as given by 

2 4 6 

Go(2 : y) = sech' y = 1 - 2 |- + 16 |- - 272 |- + .. (34) 

3 5 

Gi(2 : y) = sech^ y tanhy = y - 8 |- + 136 |- - .. (35) 

2 4 

G2(2 : y) = sech^ y tanh^ y = 2 |- - 40 ^ .. (36) 



In general the n*'' column of the triangle scheme can be related to the power series for G„_i(y). 
The numbers at the nodes in the first column of this diagram are the Tangent numbers which 
are simply related to the Bernoulli numbers. 

When cr 7^ 1, but a finite number, it is simply a scale factor and the resulting functions may be 
obtained using a scaled variable instead of y. However, when a = but Xj ^ 0, the functions 
that arise in this limit have a new structure and must be studied separately. This is the case we 
examine next. 
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2.3 Study of the case a = 

We consider the case when one of {a, f3) takes the vakie 1, the other — > oo, cr — and (ja/3 — ?► 1. 
In this hmit — >■ (1 + j), and the commutator [R,L] — )■ /, the unit matrix. The generators 
corresponding to this algebra are the raising and lowering operators of the SHO (i.e), L = a , R = 
a^. Using the relation 

2 

it^^o (sech y v^) " ^ exp(^) (37) 



in ([7]- in]) leads to the result that 



exp(zy(a + flt)) = (exp(iya^')) (exp(-^)) (exp(i?/a)) (38) 

which is well known in laser studies as the Unitary operator which may be used to generate 
the coherent states of the electromagnetic field from the vacuum state. It may be established 
that 

G„ ^ ^ exp^ (39) 

d 



„ G„+i = (V^W^ G„ - V^rT2 Gn+2) ■ (40) 
oy 

The diagram corresponding to G„ has square root factors in the arms and nodes. Simpler 
diagrams result if we consider 

G„ = = ^ exp(^) or G„ = G„ = y" exp(^) (41) 



whose power series expansion is coded by the nodes of 
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\ 
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1 






1 

\ 




2 1 
/ \ 


1 

\ 




1 
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1 




3 




1 

/ 




1 3 
\ / 


1 




2 

\ 
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1 


\ 






\ 








15 






15 




1 






1 







12 



15 15 
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where the columns on the left diagram code for the power series of Gn and the columns on the 
right diagram code for Gn- For example, 

1 + - 15|^ + .... (42) 

y - 3 |- + 15 |- .... (43) 

/2T - 2! + ^44) 



Go — 


Go 


= Go = 


Gi = 


Gi 


= Gi = 


G^ _ 


G2 





2.4 Generalisation to diagrams for particle and hole states 

The analysis of the cases considered so far indicates that for general values of a and /3 the 
functions defined by (fT2|) would correspond to functions generated from the nodes of a class of 
diagrams which have a systematic pattern of weights along the diagonals going down to the left 
or to the right. If neither of {a, (3) equals 1 the 'particle' and 'hole' states no longer decouple 
and the triangular diagrams are no longer adequate to express the functions G„ that arise from 
the general algebra. However a triangle diagram reflected about the vertical axis produces a 
diamond formation and such diagrams may be used to accomodate the feature that the nodes 
to the right of the central vertical line may be used for describing the 'particle' states and the 
nodes to the left of the central spine may be used for the 'hole' states. 





A-1 


Ao 
\ 






• 




• 






Ao 


Ai 




\ 




\ 


• 


• 






^ \ 


A-1 


Ao 
\ 


Ai 




• 




• 


A_3 A_2 


A-1 


Ao 


Ai 


\ ^ 


\ 




\ 


• 


• 






\ 


A-1 


Ao 
\ 


Ai 




• 




• 




A-1 


Ao 






\ 







\ 

A2 



It is evident from this diagram that when A_i = the weights attached to the links in the first 
column to the left of the central spine all vanish and the two halves of the diagram uncouple and 
the values at the nodes on the central spine and to the right hand side of the spine are no longer 
affected by the left side of the diagram and leads to the triangular diagrams we have considered 
earlier. Similar decouplings of the diagram into two parts will happen whenever one of the AjS 
— > 0. We next consider an example of an algebra which leads to functions coded by diagrams in 
diamond formation. 



If the parameters of the algebra are chosen so that a — (3 oo, a ^ 0, ajSa 1, in this 
limit Xj — > 1 independent of j, then S vanishes because R and L now commute and in Q 
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/ = 1, .9 = 1, and (|lll) simplifies. The resulting solutions 



y 



2j+n 



j=0 



(45) 



may be identified as the Bessel functions of integer orders (Abramowitz and Stegun 1965) that 
satisfy 

^J„+i(2y) = ( J^{2y) - Jn+2{2y)) (46) 

which is the Xj — >■ 1 limit of p^ . The power series expansion of these solutions may be 
represented by the diagram 



1 1 
\ / 
4 



1 

\ 



1 1 
\ ^ 
4 



10 



10 



20 

The central spine corresponds to n = and columns to the right of the spine correspond to 
positive integer values of n and the columns to the left of the spine correspond to negative 
integer vaues of n. The upper part of the diagram is the Pascal triangle with the nodes along 
the horizontal rows providing the binomial coefficients. We now have the additional insight that 
the numbers at the nodes of the columns may be used to provide : 



M2y) = 1 - 2 1^ + 6 if - 20 ^ 



M2y) 



y 



r 

3! 



10 



r 

5! 



1 


y' 


y^ 


0!0! " 


1!1! 


^ 2!2! 


y 


y3 


y5 


0!1! 


1!2! 


^ 2!3! 



y 

3!3! 



(47) 
(48) 



In this manner the columns of the diagram, when extended, provide the elements for constructing 
the power series for Bessel functions, the n*'' column to the right of the spine codes for J„(2y) and 
the columns to the left of the spine code for (— )" J_„(2?/). Using the properties of the binomial 
coefficients it is evident that the values at the nodes in the m*^ row add to 2™"^ when added with 
the same sign and add to zero for all rows except for the first row, if added with alienating signs. 
These features together with a property of the Bessel functions, {viz.) J-n{z) = (— )" J„(2), can 
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be used to infer from the diagram that 

oo 

.h{2y) + 2 ^ J2k('2y) = 1 (49) 
Jo(2y) + 2 ^ (-)''^ J2fe(2y) = 1 - ^ + ^ - - = cos2y (50) 

oo ^3 3 (-.^ 5 

2 ^ (-)'=+! ./2fe+i(2?y) = 2y - ^ + ^ - ... = sin2y (51) 



fc=i 



which are well known sumrules for Bessel functions (Gradshteyn and Ryzhik 1965). This is an 
example of how diagrams may be used not only to capture the power series for a whole class of 
functions but also extract sumrules satisfied by the functions. 



3 Generalisation of Unitary operator 

The normal ordered form of Ui given in ([7|) is a special case of a more general operator constructed 
as the exponential of {aL + hR + cS) where (a, &, c) are complex numbers. The normal ordered 
decomposition can be performed using a generalisation of the method used to obtain the result 
given in ([7]). For the algebra governed by © it may be shown that 



U. = 



;exp(aL + bR + cS)) = (exp(6/i?)) ( ^ b') (g^) (exp(a/L)) (52) 

j 

= (exp(a/L)) (g^) (exp(6/i?)) (53) 

i 

tang q secq j— 

J = 7 1 7 , 9 = 7 ; r , q = \/ a b a - c/ a-^ . (54) 

{ q — c a tang) [ q — ca tang ) 

For the case a ^0, a = b = iy, c = 0, C/2 reduces to Ui and - reduce to © - dH). 

The factorised form of U2 may be used to obtain the rotation matrices for the angular momentum 
algebra for which {R,L,S) ^ (^J+, Jz) and a = —\. Choice of the coefficients {a,b,c) to 
be related to the cartesian components of a vector W with polar coordinates {u!,6,(j)) by the 
relations 

a ~ uj sa\d (exp — z0) , b — i\pi uj sin(? (expi(/)) , c = — 2za; cos0 (55) 
leads to the normal and anti-normal ordered forms 

U = (exp(2iW.J)) (56) 

= (exp(a+ /i e-^^)) ( |jm) (s)"^"(jm|) (exp(a_ /i e+'^)) (57) 

= (exp(zJ_ h*e^^*)) ( |jm) (s*)+""(jm| ) ( exp(z J+ e"^*)) (58) 
a/2 

h = — sin0 sinw , s = (cosw — i cos 9 sinw) (59) 
s 
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where \jrn) are eigenstates of and Jz with eigenvahies j{j + 1) and m respectively. The matrix 
elements of the Unitary operator U may be computed for any W and J. For example, if J = 1 
then using the above representations it can be established that 







-1) " 




u 




0) 








^1) . 





h*^ e^^'t' i h*s*^ e"'' 



„*2 



|1, 0) 
II, +1) 



(60) 



which provides the matrix elements of U for arbitrary rotations W of the J = 1 system. Fur- 
thermore, ii 9 = 7r/2, = 0, then 



cos^ uj 



■ 2 

sm Lu 



-V2 si 



smw cos a; 



sin^ UJ 



\f2 sin ijj cos 

1 — 2 sin^ w \f2 sin w cos w 
— -\/2 sin a; cos w cos^ 



(61) 



4 Conclusions 



In this study we have shown that a variety of solutions emerge from the expectation values 
of a Unitary operator constructed from the generators of a SU(1,1) algebra. Coherent states 
and squeezed states of the SHO familiar in laser studies emerge as special cases. The rotation 
matrices associated with the angular momentum algebra also emerges as a special case. We have 
shown that the functions which arise as expectation values of the Unitary operartor satisfy a 
simple two term recursion relation and have power series expansions which can be captured by 
diagrams similar to the Pascal triangle familiar in the study of binomial coefficients. 

We have studied the (nO) matrix elements of the Unitary operator C/i in ([7|) in depth. For the 
general case we note that if we define 



G„„(a,/3;a;2/) = (-^)"-'" (n|C/i|m) 



(62) 



it may be shown that Gnm may be evaluated by the same steps as for G„ by starting from (|5]) 
and the resulting hypergeometric function shows that Gnm can be related to G„ in - (fT51) 
through the mapping given by 



G„jn{a,l3;a;y) = G(„_„i)o(a + to, /3 + to; ct; y) = Gn-m{a + m, {3 + m; a;y) 



(63) 



Diagrams that code for the power series of Gnm may be drawn by adopting the procedure used for 
G„ by starting in the column corrsponding to to and considering all paths proceeding downward 
with appropraite weights Xj attached to the arms. 
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